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Within an isospin- and momentum-dependent transport model using as an input nucleon density 
profiles from Hartree-Fock calculations based on a modified Skyrme-like (MSL) model, we study 
how to eliminate the influence of neutron-skin size difference of initial colliding nuclei in probing 
the nuclear symmetry energy. Within the current experimental uncertainty range of neutron-skin 
size of 208 Pb, the Pb+Pb collisions are performed in semicentral and peripheral collisions with 
impact parameters of 5 and 9fm and at beam energies from 50 MeV/nucleon to 1000 MeV/nucleon, 
respectively. It is shown that combination of neutron and proton collective flows, i.e., neutron- 
proton differential elliptic flow, neutron-proton elliptic flow difference, neutron-proton differential 
transverse flow and neutron-proton transverse flow difference, can effectively eliminate the effects 
of neutron-skin size difference and thus can be as useful sensitive observables in probing nuclear 
matter symmetry energy in heavy-ion collisions. Moreover, the combined neutron-proton stopping 
power including the neutron-proton differential stopping power and neutron-proton stopping power 
difference can also eliminate the effects of neutron-skin size difference and shows some sensitivities 
to symmetry energy especially at the lower beam energy. 

PACS numbers: 25.70.-z, 24.10.Lx, 21.65.-f 


I. INTRODUCTION 

The density dependence of nuclear symmetry energy 
E sy m{p) as one of important issue in isospin physics has 
been studied for the past few decades because of its very 
importance not only in understanding the structure of ra¬ 
diative nuclei in nuclear physicsIh-Q but also its crucial 
roles in nuclear astrophysics [^-Q- Up to now, although 
many useful experimental observables have been 

proposed to determine the nuclear symmetry energy in 
heavy-ion collisions, they have had limited success be¬ 
cause of its sensitivity not only to nuclear symmetry en¬ 
ergy but also other poorly known physical effects [20l. I2l| . 
Nevertheless, a key step in determining the nuclear sym¬ 
metry energy is the determination of experimental ob¬ 
servables which can be serve as clean and sensitive probes 
mm. Therefore, the search of the experimental ob¬ 
servables sensitive only to nuclear symmetry energy than 
others is a crucial task in determination of nuclear sym¬ 
metry energy. 

It is well known that heavy-ion reactions induced by 
neutron-rich nuclei provide an important opportunity to 
constrain the symmetry energy in a broad density range 
m~m]- To initialize transport models of heavy-ion re¬ 
actions, it is necessary to know the nucleon density pro¬ 
files for the two colliding nuclei. Generally speaking, one 
should use the Thomas-Fermi method to extract the den¬ 
sity profiles for the colliding nuclei. Practically, one usu¬ 
ally use other methods to approximate Thomas-Fermi 
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method during initializing the colliding nuclei such as 
the droplet model used in the isospin-dependent quan¬ 
tum molecular dynamics (IQMD) model [HI [2f| and 
the Skyrme Hartree-Fock model used in the isospin- 
dependent Boltzmann-Uehling-Uhlenbeck (IBUU) model 
[3d ]. In any case, the initialization should guarantee 
the numerical stabilization of colliding nuclei and fit 
the bulk nuclei properties as good as possible; see, e.g., 
Refs. j3ll. j32|. The main problem of these methods is 
the numerical stability of initial colliding nuclei in the 
subsequent reaction. However, this allows one to easily 
separate effects on the final observables due to the ini¬ 
tial state from those due to the reactions 133} . Based 
on this consideration, by examining the relative effects 
of neutron-skin size in initial nuclei and the symmetry 
energy at suprasaturation densities reached in heavy-ion 
collisions on the charged pion ratio in the final state, we 
have shown recently that the 7r _ /7r + ratio is sensitive not 
only to nuclear symmetry energy but also the neutron- 
skin size of initial colliding nuclei especially in peripheral 
collisions, see Ref. t 33|. Nevertheless, as mentioned above 
a crucial task in determination of symmetry energy is to 
find the clean experimental observables. Under this con¬ 
sideration in mind, we are going to find which experimen¬ 
tal observables are sensitive to nuclear symmetry energy 
than neutron-skin size difference of initial colliding nu¬ 
clei. To this end, in this work we investigate how to 
eliminate the influence of the neutron-skin size difference 
of initial colliding nuclei in probing the nuclear symmetry 
energy in Pb+Pb heavy-ion collisions. It can be found 
later that the combined neutron-proton collective flow 
and stopping power can effectively eliminate the effects 
of neutron-skin size difference of initial colliding nuclei 
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but keep the effects of symmetry energy especially at the 
lower beam energy. 


II. THE MODEL 

In this part we briefly describe the model used in 
the present study, i.e., the isospin- and momentum- 
dependent Boltzmann-Uehling-Uhlenbeck transport 
model [10] of version IBUU11 f+ij]. The momentum 
dependence of both the isoscalar 00 and isovector 
[30l HOHUl parts of the nuclear interaction is impor¬ 
tant in understanding not only many phenomena in 
intermediate-energy heavy-ion collisions but also ther¬ 
modynamical properties of isospin-asymmetric nuclear 
matter [4314451 ] . The mean-field potential for a nucleon 
with momentum p and isospin t can be written as [13] 
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FIG. 1: (Color online) The density dependence of the sym¬ 
metry energy. Taken from Ref. 1331 . 
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In the above, p = p n + p v is the nucleon number density 
and S = (p n — Pp)/p is the isospin asymmetry of the nu¬ 
clear medium; p n t p ) denotes the neutron (proton) density, 
isospin r is 1/2 for neutrons and —1/2 for protons, and 
f{p) is the local phase space distribution function. The 
expressions and values of the parameters A u (x). Ai(x\ 
a, B : C T<T , C T -r , and A can be found in Refs. fdol. l+il . 
and they lead to the binding energy of —16 MeV, incom¬ 
pressibility 212 MeV for symmetric nuclear matter, and 
symmetry energy E sym (p 0 ) = 30.5 MeV at saturation 
density po = 0.16 fm -3 , respectively. 

The variable x is introduced to mimic different forms 
of the symmetry energy predicted by various many-body 
theories without changing any properties of symmetric 
nuclear matter and the value of E sym (po). At suprasat- 
uration densities, although the IBUU calculations favour 
by comparing with the FOPI data a super soft symme¬ 
try energy with x = 1 (47} , compared to the FOPI data 
the ImIQMD calculations by Feng et al [48[ show a su¬ 
per hard symmetry energy. Therefore, to evaluate the 
relative effects of symmetry energy we use two values of 
x = 1 and x = 0 as the so-called soft and stiff symmetry 
energy parameters as shown in Fig. □ 0. It should be 
mentioned that the current uncertain range of symmetry 
energy at suprasaturation densities is much larger than 
the one used here 0, 0. The density dependence of 
E syrn (p ) around po is generally characterized by the slope 
parameter L = 3po(dE sym /dp) p - Po . The softer (stiffer) 
E sy m{p ) with x = 1 {x = 0) has a value of L = 16.4 
(62.1) MeV. 
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FIG. 2: (Color online) The neutron and proton density pro¬ 
files for 208 Pb with neutron-skin thickness of 0.1 and 0.3 fm, 
respectively. Taken from Ref. [33| . 


To examine effects of the neutron-skin size of initial col¬ 
liding nuclei, we initialize nucleons in phase space using 
neutron and proton density profiles predicted by Hartree- 
Fock calculations based on the MSL model [50, |5lj]. Dif¬ 
ferent values of neutron-skin thickness can be obtained 
by changing only the value of L in the MSL0 force [HU 
while keeping all the other macroscopic quantities the 
same. Shown in Fig. [2] are the density profiles corre¬ 
sponding to a neutron-skin thickness S of 0.1 and 0.3 fm 
of 208 Pb 0], which are in the range of about 0.11 ± 0.06 
fm from 7r + -Pb scattering 0 ] to 0.33lg’[g fm from the 
PREX-1 experiments using parity violating e-Pb scat¬ 
tering 0 ]. Although these available data suffer from 
large uncertainties, it was shown very recently within a 
relativistic mean-field model 0] that a neutron-skin for 
208 Pb as thick as 0.33 + 0.16 fm reported by the PREX-I 
experiment 0 ] can not be ruled out although most other 
studies have reported much smaller average values albeit 
largely overlapping with the PREX-I result within error 
bars. It is expected that the proton distributions are al- 
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FIG. 3: (Color online) Evolution of the rms radius differ¬ 
ence between neutron and proton initially setting as O.lfm 
in 208 Pb+ 208 Pb collisions with impact parameters of 5, 9, 15 
and 20fm at beam energy of 400MeV/nucleon, respectively. 


III. RESULTS AND DISCUSSIONS 


We now present results of the study in the following. 
Considering that there is no obvious neutron-skin effects 
in head-on heavy-ion reactions as shown in our previous 
work [33|, we thus carry out this study in semicentral and 
peripheral 208 Pb+ 208 Pb collisions with impact parame¬ 
ters of 5 and 9fm. On the other hand, to distinguish 
the difference of neutron-skin size of S'=0.1 and 0.3fm, 
we have performed large scale calculations with 4 x 10 5 
events in each case reported here. Thus, the statistical 
error bars are smaller than the plotting symbols in most 
plots. In addition, the percentage used in some figure is 
to avoid too much numerical digit after the decimal. 


A. The elliptic flow 


most identical, while the neutrons distribute differently 
in the two cases considered. 

Before showing the results of the studies, one must 
check the stability of the ground state nucleus. To this 
end, one can check the time evolution of the rms radius 
difference between neutrons and protons in the projec¬ 
tile and/or target with impact parameter to be infinite 
in the reaction model. This is because when the impact 
parameter is infinite, the projectile and target can not 
touch with each other and the corresponding interaction 
between them becomes zero, they are in their ground 
states and move along their initial trajectory; it is natu¬ 
rally that the initial rms radius difference between neu¬ 
trons and protons remains unchanged. Shown in Fig. [3] 
is the evolution of the rms radius difference between neu¬ 
trons and protons initially setting as O.lfm both in the 
projectile and target for different impact parameters at 
beam energies of 400MeV/nucleon, respectively. Firstly, 
it can be found that the value of the rms radius dif¬ 
ference between neutrons and protons decreases rapidly 
with the impact parameter decreasing due to the inter¬ 
actions between the projectile and target as well as the 
corresponding collisions among nucleons increasing. Sec¬ 
ondly, it is expected that the initial rms radius difference 
between neutrons and protons in the projectile and/or 
target approximately approaches stable as the impact pa¬ 
rameter increasing to be 20fm. Certainly, a long-period 
and small-amplitude oscillation of the rms radius differ¬ 
ence is still seen due to the distance between the projec¬ 
tile and target inadequate far and the tiny fluctuation 
from collisions between nucleons in the same nucleus. 
However, this level of stability of the ground state nu¬ 
cleus more or less lasts long enough, which should be 
reflected in the final reaction production. 


The elliptic flow has been widely used to study the 
properties of the hot and dense matter formed in the early 
stage of heavy-ion collisions at relativistic and interme¬ 
diate energies, see, e.g., Refs. (Hi, [H, [56j. To probe the 
density dependence of symmetry energy, we show in Fig. 
[Hand Fig. [5]transverse momentum dependence of elliptic 
flow of the midrapidity ( \y/ybeam\ < 0.5) neutrons and 
protons in semicentral and peripheral 208 Pb+ 208 Pb colli¬ 
sions with impact parameters of 5 and 9fm and at beam 
energies from 50 to 1000 MeV/nucleon, respectively. It 
can be found that the elliptic flow of nucleons shows a 
transition from in-plane to out-of-plane with the beam 
energy increasing from 50 to lOOOMeV/nucleon. This is 
because at the lower beam energies, the mean field dom¬ 
inates the reaction dynamics and causes the in-plane en¬ 
hancement of emitted reaction products, while with the 
beam energies increasing, the mean field becomes less 
important and the collective expansion process based on 
nucleon-nucleon scattering starts to be predominant, and 
the squeezed out elliptic flow as the result of shadowing 
of spectator starts to become negative (57l - [59j . However, 
the symmetry energy effects on elliptic flow of nucleons 
are not obvious as reported in previous Refs. H3, 
and the neutron-skin effects on elliptic flow of nucleons 
are comparable with or larger than that of symmetry en¬ 
ergy. Nevertheless, the elliptic flow of both neutron and 
proton has larger value with thinner neutron-skin due to 
more nucleons in spectator generating the stronger shad¬ 
owing effects on in-plane emitted nucleons and leading 
squeezed out elliptic flow of out-of-plane nucleons to be 
larger. In other words, the effects of neutron-skin size 
on elliptic flow of neutron and proton are approximately 
identical except for the additional Coulomb repulsion be¬ 
tween protons. This naturally leads us to check whether 
combining the neutron and proton elliptic flows can elim¬ 
inate the effects of neutron-skin size difference but keep 
the effects of symmetry energy. To this end, we formulate 
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FIG. 4: (Color online) The transverse momentum dependence 
of elliptic flow of the midrapidity (\y/ybeam\ < 0.5) neutrons 
in semicentral and peripheral 208 Pb+ 208 Pb collisions with im¬ 
pact parameter of 5 and 9fm and at beam energies from 50 
to 1000 MeV/nucleon, respectively. 


the so-called neutron-proton differential elliptic flow 
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where N(u ), N n (u) and N p (u) are the numbers of nucle¬ 
ons, neutrons and protons at parameter u which denotes 
the rapidity y or transverse momentum p t , and Wi is 1 
for neutrons and —1 for protons, respectively. Shown 
in Fig. [G] is the rapidity dependence of neutron-proton 
differential elliptic flow in semicentral and peripheral 
208 Pb+ 208 Pb collisions with impact parameters of 5 and 
9fm and at beam energies from 50 to 1000 MeV/nucleon, 
respectively. It is expected that the neutron-proton dif¬ 
ferential elliptic flow shows obvious sensitivity to the 
symmetry energy. This can be understandable since the 
combined elliptic flow constructively maximizes the ef¬ 
fects of the symmetry potential but minimizes the ef¬ 
fects of the isoscalar potential similar to that of neutron- 
proton differential transverse flow proposed by Li Jilj. 
On the other hand, noticing that the neutron-proton dif¬ 
ferential elliptic flow in nridrapidity ( \y/ybeam\ < 0.5) is 
hardly affected by the neutron-skin size difference but 
sensitive to symmetry energy, the transverse momentum 
dependence of neutron-proton differential elliptic flow of 
the midrapidity ( \y/ybeam\ < 0.5) is shown in Fig. [7] It 
can be found that the transverse momentum dependence 


FIG. 5: (Color online) Same as Fig. [4] but for protons. 



FIG. 6: (Color online) The rapidity dependence of neutron- 
proton differential elliptic flow in 208 Pb+ 208 Pb collisions with 
impact parameters of 5 and 9fm and at beam energies from 
50 to 1000 MeV/nucleon, respectively. 


of neutron-proton differential elliptic flow is indeed more 
sensitive to symmetry energy but hardly sensitive to the 
neutron-skin size difference especially at the lower beam 
energy. 

Another combination of the neutron and proton elliptic 
flows is the direct difference of them proposed in Ref. 
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FIG. 7: (Color online) The transverse momentum dependence 
of neutron-proton differential elliptic flow of the midrapidity 
( \y/Vbeam\ < 0.5) in 208 Pb+ 208 Pb collisions with impact pa¬ 
rameters of 5 and 9fm and beam energies from 50 to 1000 
MeV/nucleon, respectively. 


FIG. 8: (Color online) The transverse momentum dependence 
of neutron-proton elliptic flow difference of the midrapidity 
(\y/ybeam\ < 0.5) in 208 Pb+ 208 Pb collisions with impact pa¬ 
rameter of 5 and 9fm and beam energies from 50 to 1000 
MeV/nucleon, respectively. 


0,H3 defined as 


v^iPt) =< t£(pt) > - < v p 2 (p t ) >, (3) 


which should also be sensitive to symmetry energy since 
it is the special case of the neutron-proton differential 
elliptic flow, i.e., when the neutron and proton have 
the same multiplicities but different average elliptic flow. 
Shown in Fig. [8]is the transverse momentum dependence 
of neutron-proton elliptic flow difference of the midra¬ 
pidity (\y/ybeam\ < 0.5) in semicentral and peripheral 
208 Pb+ 208 Pb collisions with impact parameters of 5 and 
9fm and at beam energies from 50 to 1000 MeV/nucleon, 
respectively. It is seen that the transverse momentum 
dependence of neutron-proton elliptic flow difference is 
indeed sensitive symmetry energy but hardly sensitive 
to the neutron-skin size difference especially at the lower 
beam energy. These indicate that combination of neutron 
and proton elliptic flows, i.e., neutron-proton differential 
elliptic flow and neutron-proton elliptic flow difference, 
can effectively eliminate the effects of neutron-skin size 
difference especially at the lower beam energy and thus 
can be as useful sensitive observables in probing nuclear 
matter symmetry energy in heavy-ion collisions. 


B. The transverse flow 


The neutron-proton differential transverse flow as a 
good tracer of the symmetry potential defined as 
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was proposed to measure symmetry energy in heavy-ion 
collisions due to its advantages of combining construc¬ 
tively effects of the symmetry potential on the isospin 
fractionation and the collective flow [(31] . When the neu¬ 
tron and proton have the same multiplicities but differ¬ 
ent average transverse flow, the neutron-proton differen¬ 
tial transverse flow naturally becomes the direct differ¬ 
ence of the neutron and proton transverse flows, i.e., the 
neutron-proton transverse flow difference, 


Px p {y) =< Pxiv) > - < p p (y) > ■ ( 5 ) 

Similar to the neutron-proton differential elliptic flow and 
neutron-proton elliptic flow difference, combination of 
the neutron and proton transverse flows, i.e., neutron- 
proton differential transverse flow and neutron-proton 
transverse flow difference, can effectively eliminate the 
effects of neutron-skin size difference but keep the effects 
of symmetry energy especially at the lower beam energy 
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as shown in Fig. [9] At the beam energies of 200 and 
lOOOMeV/nucleon, the neutron-proton differential trans¬ 
verse flow and/or neutron-proton transverse flow differ¬ 
ence are less sensitive to nuclear symmetry energy com¬ 
pared to those at 50MeV/nucleon. 



FIG. 9: (Color online) The neutron-proton differential trans¬ 
verse flow (left panel) and neutron-proton transverse flow dif¬ 
ference (right panel) in 208 Pb+ 208 Pb collisions with impact 
parameter of 5 and 9fm and beam energy of 50MeV/nucleon, 
respectively. 



Beam energy (MeV) 

FIG. 10: (Color online) The beam energy and impact pa¬ 
rameter dependence of stopping power of all nucleons in 
208 Pb+ 208 Pb collisions with the stiff (upper panel) and soft 
(lower panel) symmetry energy, respectively. 


C. The stopping power 


Collective flow, generated by pressure gradient of dense 
nuclear matter formed in heavy ion collision, is closely re¬ 
lated to nuclear stopping power. Large stopping power 
can lead to a remarkable pressure gradient in the com¬ 
pressed nuclear matter. It is also believed that the stop¬ 
ping power governs most of the dissipated energy, the 
amplitude of large collective motion, the maximum at¬ 
tainable baryon and energy densities, as well as the ther- 
malization of the collision system I62l Ki ll . Two kinds 
of ratio of transverse to longitudinal quantities are com¬ 
monly used to measure the degree of stopping; one is the 
energy-based isotropy ratio, another is the momentum- 
based isotropy ratio. They are actually the different 
forms of physics realization of the classical Maxwell dis¬ 
tribution assumption (fi5i |. Here, the momentum-based 
isotropy ratio is employed to measure the degree of stop¬ 
ping, its definition is 


n _ 2 X) \Pti\ 

~~ ttE W 


( 6 ) 


wher epti(p z i) is the transverse (longitudinal) momentum 
of nucleon in center of mass system, and the sum runs 
over all products event by event. It is believed that full 
stop ping is reached when the ratio R reaches the value 
of 1 [66j; and the superstopping of the ratio R > 1 is 
explained by the preponderance of momentum flow per¬ 
pendicular to the beam direction [67| . 
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FIG. 11: (Color online) The rapidity dependence of stopping 
power of neutrons in 208 Pb+ 208 Pb collisions with impact pa¬ 
rameter of 5 and 9fm and beam energies from 50 to 1000 
MeV/nucleon, respectively. 
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FIG. 12: (Color online) Same as Fig. QT]but for protons. 
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FIG. 14: (Color online) The transverse momentum depen¬ 
dence of neutron-proton differential stopping power (left 
panel) and neutron-proton stopping power difference (right 
panel) of the midrapidity ( \y/ybeam\ < 0.5) in 208 Pb+ 208 Pb 
collisions with impact parameter of 5fm and beam energies 
from 50 to lOOOMeV/nucleon, respectively. 


energy, respectively. It is seen the stopping power of all 
nucleons is decreasing with the beam energy and impact 
parameter increasing regardless of the stiff or soft symme¬ 
try e nergy . These are consistent with previous results in 
Ref. [63ll64| . However, the stopping power of all nucleons 
is not sensitive to symmetry energy at all. Is the stop¬ 
ping power of neutrons or protons sensitive to symmetry 
energy? Can the combined stopping power eliminate the 
influence of neutron-skin size difference of initial colliding 
nuclei but keep the effect of symmetry energy? To an¬ 
swer these questions, let’s first define the neutron-proton 
differential stopping power R np and neutron-proton stop¬ 
ping power difference R n ~ p similar to those of combined 
collective flow as follows, 


R np (u ) 


N n (u) 

N(u) 


< R n (u) > 


N p (u) 

N(u) 


< R v (u) >, 


(7) 


FIG. 13: (Color online) The rapidity dependence of neutron- 
proton differential stopping power in 208 Pb+ 208 Pb collisions 
with impact parameter of 5 and 9fm and beam energies from 
50 to 1000 MeV/nucleon, respectively. 


Let’s first look at the global dependence of the stop¬ 
ping power on beam energy, impact parameter and sym¬ 
metry energy for a given neutron-skin thickness of O.lfrn. 
Shown in Fig. [TO] is the beam energy and impact pa¬ 
rameter dependence of stopping power of all nucleons in 
208 Pb+ 208 Pb collisions with the stiff and soft symmetry 


R n ~ p (u) =< R n {u) >^<R p (u)>, (8) 

where N(u ), N n (u) and N p (u) are the numbers of nucle¬ 
ons, neutrons and protons at parameter u which denotes 
the rapidity y or transverse momentum p t , respectively. 
Shown in Figs. fTTl fl2l and fT3l are the rapidity dependence 
of stopping power of neutrons, protons and neutron- 
proton differential stopping power in 208 Pb+ 208 Pb colli¬ 
sions with impact parameters of 5 and 9fm and at beam 
energies from 50 to lOOOMeV/nucleon, respectively. It 
can be found that although the stopping power of neu¬ 
trons and/or protons does not show obvious sensitivities 
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FIG. 15: (Color online) Same as Fig. [14] but with impact 
parameter of 9fm. 

to symmetry energy, the neutron-proton differential stop¬ 
ping power shows obvious sensitivities to symmetry en¬ 
ergy especially for midrapidity nucleons and at the lower 
beam energy. In addition, the neutron-proton stopping 
power is hardly affected by the neutron-skin size differ¬ 
ence of initial colliding nuclei. This is very similar to that 
of combined transverse and elliptic flows due to their cor¬ 
relation as pointed out in Ref. [621 ]. i.e., flow is generated 
by pressure gradients established in compressed matter, 
while the achieved density is connected to the degree of 
stopping. 

Finally, noticing that the midrapidity neutron-proton 
differential stopping power is obvious sensitive to sym¬ 
metry energy but hardly affected by the neutron-skin 
size difference of initial colliding nuclei, the transverse 
momentum dependence of midrapidity neutron-proton 
combined stopping power including the neutron-proton 
differential stopping power and neutron-proton stopping 


power difference is shown in Figs. [l4]and[T5l respectively. 
It is shown that the transverse momentum dependence 
of midrapidity neutron-proton combined stopping power 
shows more obvious sensitivities to symmetry energy but 
hardly affected by the neutron-skin size difference espe¬ 
cially at the lower beam energy. 

IV. SUMMARY 

In this work by studying the influence of the neutron- 
skin size of initial colliding nuclei on the collective 
flow and nuclear stopping power, we have showed 
how to eliminate the influence of the neutron-skin size 
difference of initial colliding nuclei in semicentral and 
peripheral Pb+Pb collisions at beam energies from 50 
to 1000 MeV/nucleon, respectively. Noticing that the 
effects of neutron-skin size on collective flow of neutron 
and proton are approximately identical except for the 
additional Coulomb repulsion between protons, it is 
thus that combination of neutron and proton collective 
flows, i.e., neutron-proton differential transverse and 
elliptic flows and neutron-proton transverse and elliptic 
flow differences, can effectively eliminate the effects of 
neutron-skin size difference especially at the lower beam 
energy and thus can be as useful sensitive observables of 
nuclear matter symmetry energy in heavy-ion collisions. 
In addition, the combined stopping power including 
neutron-proton differential stopping power and neutron- 
proton stopping power difference also shows some 
sensitivities to symmetry energy but hardly affected 
by neutron-skin size difference of initial colliding nuclei 
especially at the lower beam energy. 
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